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ABSTRACT 

We investigate the radio and 7-ray beaming properties of normal and mil- 
lisecond pulsars by selecting two samples from the known populations. The first, 
Sample G, contains pulsars which are detectable in blind searches of 7-ray data 
from the Fermi Large Area Telescope. The second, Sample R, contains pulsars 
detectable in blind radio searches which have spin-down luminosities E > 10^^ erg 
s~^. We analyse the fraction of the 7-ray-selected Sample G which have de- 
tectable radio pulses and the fraction of the radio-selected Sample R which have 
detectable 7-ray pulses. Twenty of our 35 Sample G pulsars have already ob- 
served radio pulses. This rules out low-altitude polar-cap beaming models if, 
as is currently believed, 7-ray beams are generated in the outer magnetosphere 
and are very wide. We further find that, for the highest-i? pulsars, the radio 
and 7-ray beams have comparable beaming factors, i.e., the beams cover similar 
regions of the sky as the star rotates. For lower- ii^ 7-ray emitting pulsars, the 
radio beams have about half of the 7-ray sky coverage. These results suggest 
that, for high-i? young and millisecond pulsars, the radio emission originates in 
wide beams from regions high in the pulsar magnetosphere, probably close to the 
null-charge surface and to the 7-ray emitting regions. Furthermore, it suggests 
that for these high-ii^ pulsars, as in the 7-ray case, features in the radio profile 
represent caustics in the emission beam pattern. 

Subject headings: pulsars: general 
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Introduction 



The observed pulses from pulsars are interpreted as an emission beam, fixed relative to 
the rotating neutron star, sweeping across the Earth. The magnetospheric emission-region 
structure, the viewing angle with respect to the rotation axis and the angle between 
the magnetic and rotation axes define the shape of the observed pulse profile. Pulsar 
radio emission is generally thought to o riginate within a few perc ent of the light-cylinder 



radius (-Rlc) above the polar cap s (e.g. 



a "radius-to-frequency" mapping (ICordes 



Lyne fc Manchester! il988i ). Many pulsars exhibit 



19781 ) which is indicative of polar cap emission. 



The increasing pulse width at lower radio frequencies is interpreted as emission at 
increasing heights above the poles along more widely directed field lines. The radio 
profiles of many young pulsars and millisecond pulsars (MSPs) have an "interpulse" , 
a second pulse component separated from the main pulse by close to 180° of pulse 
phase, or widely separated pulse components. In most cases, such pulsars do not exhibit 
radius-to- fre quency mapping, that is , they have frequency-independent pulse component 



separations fJHankins fc Fowler 



19861 ) 



The known population of 7-ray pulsars has recently been grea tly expanded by th e 



2010bl . 



six-month data release of the Fermi Large Area Telescope (LAT) (jAbdo et al. 

hereafter AlO). This sample of 46 7-ray pulsars includes both young pulsars and MSPs, all 

with high spin-down luminosities: 



E 



p3 



> 10^^ ergs \ 



(1) 



where / = 10^^ gm cm ^ is the assumed neutron star moment of inertia, P is the rotation 
period in seconds and P is the period time- derivative. 



Recent studies indicate that 7-ray emission from pulsars originates at hig! 
in the pulsar magnetosphere, near i?Lc aiid from the last open field lines (e.g.. 



altitudes 



Abdo et al 
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2009d ]bl. AlO). Two classes of high-altitude 



gener ally considered: slot-gap accelerators ( iMuslimov fc Harding 



particle acceleration and e r nission model ar e 



2004 : 



Dvks fc Rudak 



20031) from the polar cap outwards along the last ope n field lines, and outer-gap accelerators 



( jRomani fc Yadigaroglu 



1995 



Cheng fc Zhang 



19981 ) where the last closed field lines cross 



the null-charge surface defined by fi • B = (il is the rotation velocity vector, and B is 
the magnetic field vector). These models result in what we shall term wide-beam emission 
geometries. Low-altitude models for 7-ray emission have been largely ruled out by the 
generally wide and multi-component 7-ray pulse profiles and the existence of pulsed 7-ray 
emission at energies greater than a few GeV; at lo w altitudes such emi s sion i s precluded 
by electron-positron pair absorption. Furthermore, iRomani fc Watterd (120091 ) find that 



outer-gap models are statistically preferred over slot-gap models based on fits to LAT pulse 
profiles for a sample of bright 7-ray pulsars. 

Four of the five brightest 7-ray pulsars detected b y the Energetic Gamma-Ray 



Experiment Telescope (EGRET) are also radio pulsars (JThompson 



20081). This, together 



wit h the simila r radio and 7-ray pulse morphology for many young pulsars and MSPs, 



led 



Manchesten (120051 ) to suggest that, for these high-i? pulsars, the radio emission is also 



emitted in wide beam s from the outer magnetosph ere n ear -Rt;C- Within the framework 



of polar cap emission. 



Johnston fc Weisbergj ( l2006l ) and 



Karastergiou fc JohnstonI (12007 



hereafter KJ07) proposed that the radio emission from young pulsars is emitted at 
intermediate altitudes of ~ 1000 km or ~ 0.2i?LC in order to account for wide profiles of 
young pulsars. 

In this letter, we identify samples of 7-ray-selected and radio-selected pulsars and 
consider the pulsars in each sample that also emit radio and 7-ray pulses respectively. We 
then discuss the implications of our analysis for the radio beams of young pulsars and 
MSPs. 
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2. The 7-ray and radio pulsar samples 

We examine two pulsar samples: a 7-ray-selected sample (Sample G) and a radio- 
selected sample (Sample R). Each sample consists of a set of pulsars detectable using "blind 
search" techniques in the relevant band. 



2.1. The 7-ray-selected Sample G 

AlO presented a catalogue of 7-ray pulsars detected using two different techniques. 
First, 29 7-ray pulsars were found using timing epheme rides, mostly froni contemporaneous 
radio observations of 218 pulsars with E > lO^^ergs"^ (jSmith et al.ll2008l ). Twenty of these 
were young pulsars, and nine were MSPs. The second method involved a blind search for 
7-ray pulsars at the positions of ~300 unidentified 7-ray sources in the six-month LAT 
all-sky image. This search was sensitive to periodicities between 0.15 s and 2 s, and yielded 
16 7-ray pulsar discoveries. A further two pulsars with existing 7-ray rotation ephemerides 
- Geminga and PSR J1 124— 5916 - were also re-detected. While Geminga is well known for 



its lack of radio pulses ( jKassim fc Lazio 



19991), PSR J1124-5916 is a radio pulsar for which 



a 7-ray rotation ephe meris was used bec ause the extremely weak radio emission makes 



radio timing difficult ( jCamilo et al 



20021) 



Sample G consists of all 7-ray pulsars found in the LAT blind search. We further 
include those pulsars found using radio ephemerides that have phase-averaged 7-ray fluxes 
above the blind search sensitivity limit |j This is a conservative estimate, as some of the 
LAT blind search pulsars are in fact below the limit. Our inclusion in Sample G of pulsars 
where the 7-ray pulses were flrst detected using radio-based ephemerides is justifled because 



^The LAT blind search sensitivity limit is plotted as a function of the diffuse 7-ray 
background in Figure 10 of AlO. 
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the blind search over the unidentified 7-ray point sources was carried out after the search 
using radio ephemerides. In total, we included the 18 pulsars found in the 7-ray blind 
search and 17 pulsars found using radio ephemerid es, of which th r ee are MSPs. Searches for 



radio pulses from 15 of the blind search pulsars by 



Camilo et al. 



( l2009bl ) have revealed two 



detections. We do not discard t h e three blind search pulsars that have not been searched for 



radio emission by 



Camilo et al. 



(l2009bl ) because their positions have already been searched 



19851 ) with a limiting fiux 



by the Greenbank low-luminosity pulsar survey (iDewey et al. 

density of 2 niJy at 400 MHz . Therefore, of the 35 7-ray pulsars in Sample G, 20 are radio 

pulsars. 



2.2. The radio-selected Sample R 



Sample R was formed from the 218 high-E' radio pulsars observed by the LAT, of 



which 25 are also 7-ray pulsars. E values were corrected for the Shklovskii effect (IShklovskii 



I970I ) using values of P, proper m otion and NE200 



in the ATNF Pulsar CatalogucI ( jManchester et al. 



distances (ICordes fc Lazio 



20021 ) given 



20051 ): only pulsars with corrected 



E > 10'^^ ergs ^ were included in Sample R. We note that the Sample G pulsars have 
already had the Shklovskii corrections applied (AlO). 

PSR B 1509 — 58 was not detected in the LAT data, but was found by the AGILE 



7-ray telescope (jPellizzoni et al. 



2OO9I ) and so we include it in Sample R. Another pulsar , 



2010ah . 



PSR J0034— 0534, was detected using 13 months of Fermi LAT data (lAbdo et al.l 
Fourteen of the high-i? pulsars were not detected i n blind radio survey s, but were found in 



deep searc hes of supernova ren mants (SNRs) (e.g. 



Camilo et al. 



boxes (e.g. 



O'Brien et al. 



2009aj) and EGRET error 



20081 ). We discard these 14 pulsars from Sample R. This leaves 



^http : / /www. atnf . csiro.au/research/pulsar/psrcat 
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200 radio-selected pulsars in Sample R, of which 17 are 7-ray pulsars. Table [1] summarizes 
the two samples. 

3. Implications for the radio beams 

We consider the radio and 7-ray beams as being defined by emissivity functions 
Fi{6,(f)), where the index i = (r, g) indicates, respectively, the radio and 7-ray bands. The 
co-latitude angle 6 is measured from the rotation axis and (p is the longitudinal angle. The 
"latitude coverage", Ci{6), of a beam is defined as the longitude integration of Ff 

Cm= [^ F,{9,(P)d<p. (2) 

The fraction of the celestial sphere swept out by each beam - the "beaming fraction" - is 
given by 

fi = — / Fi{e,(j))smed(j)de = — Ci{e)smede (3) 

4vr Jo Jo 47r Jq 

where the normalisation is such that < /j < 1. 

To assist with visualisation, we approximate the beams by "top-hat" functions where 
Ci is constant (and non-zero) over some beamwidth Bj and zero elsewhere as illustrated 
in Figure [TJ Bj is effectively the equivalent width of the beam in latitude. We note that 



the beaming 



Watters et al 



raction /j defined by Equation [3] corresponds to the fn term defined in 



(120091 ) for such a top-hat beam. There are four possible configurations for a 
pulsar emitting in both the radio and 7-ray bands: 1. The radio and 7-ray beams have no 
overlap. 2. Partial overlap of the beams. 3. The radio beam is narrower than the 7-ray 
beam and entirely overlapping. 4. The 7-ray beam is narrower than the radio beam and 
entirely overlapping. 

If we assume that pulsar rotation axes are randomly oriented in space and that pulses 
are detectable for observer lines of sight within the beams, then the number of pulsars 



detectable in band i is proportional to the beaming fraction /j. Furthermore, provided 
the beam co-latitude di is not too small, then the number of pulsars detectable in band i 
is approximately proportional to the beamwidth Bj. Importantly, the number detected in 
both bands is proportional to the width of the overlap region Be- 

The issue of pulsar detectability is of course critical. It is well-known that the 
7-ray luminosity of a pulsar, Lg, is strongly dependent on its spin-dowii luminosity E. 



While some 7-ray emission mo dels predict Ln (x E ^ (e.g. 



empirical relation Lg oc E after iSmith et al. 



Harding 



198ll ). we assume an 



( 120081 ). This implies the commonly used 7-ray 



detectability metric of the pulsar spin-down flux at the telescope, Ed'"^, where d is the 
pulsar distance. In contrast, the radio luminosity of pulsars is essentially uncorrelated with 
both E and Ed^^. Given that nearly 2000 radio pulsars are known compared to just 46 
7-ray pulsars, it is reasonable to assume that selection effects are much more significant in 
the 7-ray band and that all 7-ray pulsars are emitting potentially detectable radio beams. 

Therefore, for the 7-ray-selected Sample G, the radio fraction, i.e., the ratio of the 
number of radio detections to the total, J-'r^c = 20/35, is an estimate of the ratio Bc/Bp. 
For configuration 3 in Figure [H Be = B^; in all other configurations Be < B^. Therefore, 

f ^ 1^ > ^r.G ^ 0.571. 

Since the values of J-'r,G follow a binomial distribution, the standard error on this mean 

value is 

a = ./^SZ^ = 0.084, 
V n 

where n = 35 is the number of pulsars in the sample. 

We conclude that, on average for 7-ray emitting pulsars, the radio beaming fraction is 
at least half the 7-ray beaming fraction. It is possible that some of the current Sample G 
pulsars will be detected in the radio band, which would increase J>,g- With the continued 
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operation of the LAT, it is certain that the sample of 7-ray-selected pulsars will increase. 
The value of J-'r^c iiiay either increase or decrease, but its uncertainty will decrease. The 
longitudinal geometries of the beams are not constrained; therefore the radio and 7-ray 
emission regions need not be co-located. 

To quantify the true latitude coverage of the radio beams with respect to the 7-ray 
beams, as opposed to just constraining the overlapping area, we compare the radio-selected 
Sample R with Sample G. First, given that the LAT is only sensitive to pulsars with high 
values of Ed~'^, the 7-ray fraction in Sample R, J^g,R = 17/200, is not significant to this 
discussion. Instead, we order the pulsars in the two samples by decreasing E and consider 
the ratios J-'j.^g and J^g,R as functions of E. For each pulsar, we compute the ratios for the 
sub-sample of pulsars with E greater than or equal to the value for that pulsar. 

The results are shown graphically in Figure [2l It is striking that the seven highest-E" 
Sample G pulsars are detected in both bands. Within the same E range, all four Sample R 
pulsars are also 7-ray emitters. It is therefore evident that J^r,G and J^g,R have comparable 
values in this high-ii^ regime. For lower- ii^ values and correspondingly lower- E'rf"^ values in 
Sample R (shown in the bottom panel of Figure [2]) Jg r reduces sharply, but this is entirely 
due the Edr"^ selection for 7-ray pulsar detections. However, it is significant that some 
radio-only pulsars with E between 2 x 10^^ and 4 x 10^^ ergs~^ have high values of Edr"^. 
These pulsars probably have 7-ray beams not directed toward us, suggesting that the radio 
beams for these high-E' pulsars are of comparable width to the 7-ray beams and do not 
always overlap with their latitude coverage. Any of configurations 1, 2, or 4 in Figure [T] 
could be applicable to this case. 

In contrast, under the assumption that radio pulse detectability is independent of E, 
Sample G is free of these selection effects and the drop in J-'r.c at E' ~ 6 x 10^^ ergs~^ can be 
interpreted as an evolution in J-'j.^g from unity for high-ii^ 7-i'ay pulsars to ~ 0.45 for lower- ii^ 
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7-ray pulsars. In general, the highest-i? pulsars are among the youngest known pulsars, as 



20001 ) ■ Furthermore, pulsars evolve 



determined, for example, by SNR associations (iKaspil 
toward lower E values. It therefore appears that, for the highest-ii^ pulsars, the radio and 
7-ray beams have comparable beaming fractions with largely overlapping radio and 7-ray 
beams. As these pulsars age, the radio beams evolve to overlap about half the latitude 
coverage of the 7-ray beams. Ultimately the pulsars become undetectable in the 7-ray band 
and the radio beams become e ven narrower as they evolve toward the beaming fractions of 



0.1 typical for older pulsars (JTauris fc Manchester 



19981). 



4. Millisecond pulsars 



J 



The lower- ii^ pulsars in Samples G and R include five MSPsJj marked with the letter 
'M' in Figure [21 Despite the differences between normal pulsars and MSPs (e.g., surface 
magnetic field strength and rotation period) we do not isolate MSPs from our samples for 
three reasons. First, the three MSPs in Sample G were 7-ray-selected, which implies some 
commonalit i es with you ng 7-ray pulsars. Se cond, the nine MSPs detected with the LAT by 



AbdoetaL 



fl2009ah and 



AbdoetaL 



(J2010al ) have values of Ed'"^ that are among those of 
the top 18 MSPs searched, which gives a 7-ray fraction of 0.5, consistent with the lower 
bound on th e radio fraction for the lower- £" Sample G pulsars (~ 0.45). Third, it has been 



suggested by IManchesterl (120051 ). based on pulse profile and giant pulse characteristics, that 
millisecond pulsars have wide radio beams similar to those of young pulsars. Therefore, in 
the discussion below on the radio beams of 7-ray emitting pulsars, we treat young pulsars 
and MSPs together. 



^While 



AbdoetaL 



( l2009al ) report the 7-ray detections of eight MSPs, three did not fall 



within either of the cutoffs applied to Samples G or R. 
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5. Discussion 



To test the consistency of our results with various empirical polar-cap models for 
pulsar beaming, we carried out a series of Monte Carlo simulations for the Sample G 
pulsars to calculate their values of J^r,G- Sample R was not considered because of the 
Ed'"^ selection effects discussed above. We assume that the 7-ray beams sweep the entire 
celest ial sphere, that is, w e assume fg = 1, consistent with outer- magnetosphere models 



[e.g. 



Watters et al. 



20091 ) ■ For each pulsar, 10^ pairs of a and ( were randomly chosen from 



sinusoidal distributions peaked at it/2 [a is the angle between the rotation and magnetic 
axes and ( is the angle between the rotation axis and the line of sight). The radio beam 
opening half-angle p, measured from the magnetic axis, was calculated for each Sample 
G pulsar and beaming model. The beaming fraction fr was then derived by finding the 
fraction of points with |a — CI < P- Since both a and ( ranged between and 71/2, emission 
from both magnetic poles was accounted for. 



We first consider the model of 



RankinI ( 19931 ) with p = 5.8P ^/^ for which we obtain 



an average beaming fraction of fr = 0.46 for the Sample G pulsars, somewhat lower than 



our estimate of fr/ fg ^ fr ^ 0.57. We no te that our assumption o : 



this co ntext; other similar simulations by 



Romani fc Yadigaroglul (119951 ) and 



fg = 1 is conservative in 



Harding et al. 



(I2OO7I ) using theoretical predictions of fg based on the outer-gap model find fr = 0.3 and 



fr = 0.1 respectively, much less than the observed value. 

We then considered the model of KJ07 which, for young pulsars, has a radio emission 
height above the polar caps of ~ 1000 km, along with beaming along the last open field 
lines. This model results in fr = 0.74, larger than our lower bound of ~ 0.57. Hence, the 
KJ07 radio beaming model is consistent with our results. However, the fact that the seven 
highest-i? pulsars in Sample G are all observed to emit radio pulses (i.e., have fr ~ 1) is 
not predicted by the KJ07 model. The product of the radio beaming fractions derived from 
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the KJ07 model of these seven pulsars is 0.22; hence, the KJ07 model has a 22% chance of 
reproducing this result. 

In contrast to the KJ07 model, outer-magnetosphere wide-beam models for the radio 
emission predict /^ ~ 1, consistent with our result that J^r,G is near unity for the highest-ii^ 
pulsars. Besides this, the re are other comp elling reasons to consider wide-beam models for 



these high-E' pulsars (cf.. 



Manchester 



20051 ) . First, there is an excess of pulsars with widely 



se parated pulse components or interpulses among young, high-i? pulsars (e.g., see Fig. 1 



of 



Weltevrede fc Johnston 



20081 ) . Many of these pulsars have pulse profiles very similar to 
observed 7-ray profiles, with bridge emission between the main components. Furthermore, 
they have frequency-independent component separations, inconsistent with low-altitude 



Narayan fc Vivekanandl (119831 ). polarization 



polar-cap models. As was first pointed out by 
results from many young pulsars indicate beams elongated in the latitude direction, with 
impact parameters, \a — C|, that are much greater than observed pulse widths. This is 
inconsistent with polar-cap models but fully consistent with outer-magnetosphere models. 

We therefore propose that the radio emission region for young high-£' pulsars is located 
high in the pulsar magnetosphere, relatively close to the null-charge surface and to the 
7-ray-emitting region. The evolution in J^-^q with E suggests a slot-gap model for the radio 
emission and hence that radio profile components represent caustics just as for the 7-ray 
profiles. This interpretation is also consistent with a transition to low-altitude polar-cap 
emission for older pulsars. 

Not all 7-ray pulsars are seen in the radio and not all high-E'c?"^ radio pulsars are seen 
in 7-rays. This implies that, while the radio and 7-ray emission regions may be in the same 
general region of the pulsar magnetosphere, they are not necessarily co-located. In the case 
of PSR J0034— 0534 however, the radio and 7-ray profiles are almost identical and align ed 



in phase, suggesting that the emission regions are in fact co-located (jAbdo et al. 



2mm- 
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6. Conclusions 

By comparing tlie relative occurrence of radio and 7-ray pulsed emission in 7-ray- 
selected and radio-selected samples, we find that the radio beaming fraction /,. appears 
to be close to unity for the highest--^ pulsars, decreasing to ~ 0.5 for lowei-E 7-ray 
pulsars. Our estimated lower bound fr, averaged over all 7-ray-selected pulsars, of ~ 0.57 
is inconsistent with low-altitude polar-cap models for the radio emission but is consistent 
with KJ07's intermediate-altitude polar-cap model. However, the KJ07 model is marginally 
inconsistent with our result that /,. ~ 1 for the highest-ii^ pulsars. With other evidence 
from radio and 7-ray pulse morphologies, these results suggest that, for high-i? pulsars, 
the radio emission originates in wide beams from the vicinity of the null-charge surface, 
possibly with a slot-gap configuration, with profile components representing caustics in the 
emission pattern. 

Improved statistics from further 7-ray and radio discoveries will help to test 
our conclusions. Modelling of radio pulse profiles with slot-gap or other wide-beam 
emission geometries for both young high-i? pulsars and MSPs would be valuable. If 
outer-magnetoshpere wide-beam radio emission from these pulsars is confirmed, it would 
have profound implications for pulsar population and evolution studies as well as our 
understanding of radio and 7-ray pulse emission mechanisms. 

We thank Simon Johnston and Kyle Watters for valuable discussions and an anonymous 
referee for constructive suggestions. GH is supported by an Australian Research Council 
QEII Fellowship (project #DP0878388). 
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Table 1. The 7-ray-selected and radio-selected samples 



Sample 7-ray pulsars Radio pulsars 



7-ray-selected (G) 35 20 

Radio-selected (R) 17 200 
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v> 



v> 



0X0r 





0c=®g<®r 



Fig. 1. — Four possible configurations for tlie radio and 7-ray beams from a pulsar emitting 
in both bands. The sectors represent the latitude coverage of the sky as the pulsar rotates for 
the 7-ray beam (G) and the radio beam (R) where the beams are approximated by "top-hat" 
functions for clarity. The beams are centered at co-latitudes 9g and 0^. and have equivalent 
widths of Qg and Qr for the 7-ray and radio bands respectively. The angular zone common 
to the two beams is shown as a hatched area and has a width ©c. If the observer's line of 
sight lies within the hatched area, both radio and 7-ray pulses can be observed. 
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Fig. 2. — The radio-detected fraction of Sample G, J^^^q (top), and the 7-ray detected fraction 
of Sample R, J^g,R (middle), computed for each pulsar for the sub-sample of all pulsars having 
E greater than or equal to the value for that pulsar. Note that E increases to the left on 
the plot. The bottom plot shows values of Ed~'^ for the Sample R pulsars. 



